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ONE-DIMENSIONAL  TWO-PHASE  CONSTANT  LAG  FLOW  IN  / 1* 

COMBUSTION  CHAMBER  OF  SOLID  PROPELLANT  ROCKET  MOTORS 


Ch'ang  Hsien-ch'i 
Abstract 


In  this  paper,  one-dimensional  two-phase  constant  lag  flow  in  combustion  chamber  of 
solid  propellant  rocket  motors  is  described.  Introducing  a  specific  heat  ratio  1C  and  a 
gas  constant  R  for  constant  lag  flow,  fundamental  relationships  analogous  to  those  of 
pure  gas  flow  are  derived.  The  effect  of  particle  velocity  lags  on  pressure-time  curves  and 
flow  field  in  chamber.  is_  analyzed^. 


Main  Notations 


Combustion  area  of  grain 
Cross-sectional  area  of  passage 
Mechanical  equivalent  of  heat 
Area  of  nozzle  throat 
Coefficient  of  combustion  rate 
Characteristic  velocity  of  propellant 
Particle  specific  heat 

Constant-pressure  specific  heat  of  gas  phase 
Constant-volume  specific  heat  of  gas  phase 
Gravitational  acceleration 
Enthalpy  per  unit  weight 

Total  enthalpy  of  1  kg  of  products  of  combustion 

Throat-passage  ratio 

Specific  heat  ratio  of  gas  phase 

Lag  constants 

Surface-throat  ratio 

Length  of  grain 


This  paper  was  received  on  November  29,  1982. 
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m  -  particle  mass  flow 
P 


m 

mc 

n 

P 

P 


r 

r. 


S 

t 

T 

TC 

v 

W 


x 

X 

p 

Pr 

PmP 


e 

P# 


Mass  flow  of  mixture  of  two  phases 

Mass  flow  of  gas  phase 

Index  of  pressure 

Pressure 

Total  pressure 

Prandtl's  number 

Combustion  rate 

Radius  of  particle 

Gas  constant  of  gas  phase 

Perimeter  of  passage 

Time 

Temperature 

Total  temperature  of  gas  phase  at  the  head  of  the  grain 

Flow  rate 

Jet  velocity 

Axial  coordinate 

Relative  axial  coordinate  i— y 

Density 

Density  of  propellant 

Density  of  Al-O, 

23  ._<»» 

Partial  particle  mass  flow,  m 

Coefficient  of  thermal  conduction  of  gas 

Coefficient  of  dynamic  viscosity  of  gas 


Subscripts 
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g  -  Gas  phase 

0  -  Cross-section  of  grain  head 

t  -  Cross-section  of  nozzle  throat 

P  -  Particle  phase 

L  -  Cross-section  of  tail  end  of  grain 


I. 


Foreword 


AJ^O^  particles  are  formed  when  composite  solid  propellants 
containing  aluminum  undergo  combustion.  Their  fractional  weight 
may  reach  30-40%.  Thus,  there  is  a  mixture  of  two  phases  (gas 
and  particle)  in  the  products  of  combustion  that  flow  in  the 
combustion  chamber  and  the  nozzle. 

Most  of  the  current  studies  on  two-phase  flows  in  solid 

propellant  rocket  engines  focus  on  the  nozzle.  Kliegel  studied 

one-dimensional  two-phase  constant  lag  nozzle  flow,  and  obtained 

[  2 1 

a  solution  similar  to  that  for  isentropic  flow  . 

The  two-phase  flow  in  the  combustion  chamber  is  one  which 
involves  mass,  and  possesses  some  new  characteristics  in 
comparison  with  that  in  the  nozzle.  In  this  paper,  we  describe 
one-dimensional  two-phase  constant  lag  flow  in  the  combustion 
chamber,  and  analyze  the  effect  of  particle  velocity  lags  on  the 
pressure-time  curves  and  the  flow  field  in  the  combustion 
chamber.  This  is  helpful  for  accurately  predicting  the  pressure¬ 
time  curve,  and  provides  more  precise  boundary  conditions  for 
computing  two-phase  nozzle  flows. 

II .  Governing  Equations 

Assume  that: ^ ^ 

1.  The  flow  is  one-dimensional  and  pseudo-steady-state. 

2.  The  friction  and  heat  conduction  between  the  flow  and  the 
passage  walls  are  neglected. 

The  Al-2^3  Particles  are  uniform  in  size,  and  are  in  fluid 
state.  The  spatial  volume  occupied  by  these  particles  as 
well  as  the  effect  of  their  Brjmnian  motion  on  pressure  are 
neglected. 


3. 


4.  The  gas  phase  is  an  ideal  gas  with  frozen  components.  It 
has  no  viscosity  except  where  it  is  in  contact  with  the 
particles . 

5.  The  force  exerted  by  mass  is  neglected. 

Accordingly,  we  have  derived  from  Ref.  [1]  the  governing 
equations  for  the  two-phase  mixture  in  the  grain  passage  of 
constant  cross-section: 

Mass  equation 

dift  » 
dx~fiTrS 


Momentum  equation 


(2) 


Energy  equation 

■^[*»  (  ft#  +'4«'fj)  +  *'(  hp +AaTg)\'m  prrSHs 


for 

From  and  ,*  ,  taking  Eq.  (1)  into  consideration, 

we  can  write  Eq.  (3)  as 

^r|[(l-e)(A,+^p^-)+e(ft,+^)-//,J-«.}-0  (4) 

Integrating  Eqs.  (2)  and  (4)  (when  x  =  0 ,  v^  =  vp  =  0 ,  /3_ 

m  =  0) ,  we  obtain,  respectively,  the  integral  form  of  the 
momentum  equation  and  the  energy  equation  of  the  two  phase 
mixture : 


Pi-P-PjO1, 


(5) 
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III.  Fundamental  Relations  of  Constant  Lag  Flows 


Define 


J9 


L (0«s£<1) 


Then,  particle  velocity  lag 


particle  temperature  lag 


Vt-Or 


t,-t. 


The  so-called  constant  lag  flow  refers  to  the  two-phase 
flow  in  which  K  and  L  remain  constant  along  the  length  of  the 
passage. 

T, 

1 .  Temperature  Ratio 

When  the  temperature  of  the  particles  exceeds  2318°K, 

the  specific  heat  of  A^O^  is  a  constant,  and 

^=»A,w^C,(7',-r,,)w. 

In  the  above  relation,  hpm  is  the  enthalpy  of  liquid  Al^O^  at  the 
temperature  Tpm  -  2318 8K. 

Denote  the  total  enthalpy  of  1  kg  of  the  gas  phase  and  that 
of  1  kg  of  the  particle  phase  by  hgg  and  hgp,  respectively.  Then 
the  total  enthalpy  of  1  kg  of  the  two-phase  mixture  is 

-  5  - 


From  Assumption  2  above  we  know  the  K_  remains  constant  along  the 
length  of  the  passage.  Let  us,  therefore,  express  Hg  in  terms  of 
the  parameters  (f  =  vp  =  0,  T^  =  Tp  =  TQ)  at  the  grain  head. 

Thus , 

#*=-(! -OCMr. +  +C, (7*. 


Inserting  the  expressions  for  Hg  and  hp  above,  as  well  as 
hg  =  CPgTg  ^nto  E<3*  >  an<^  taking  note  of  Eqs.  (7)  and  (8)  ,  we 

obtain  the  energy  equation  for  the  constant  lag  two-phase 
mixture : 


(11) 


In  the  above  equation. 


*  (1  —e)Cpf  +eCiL  r 

C'~  - C" 


1 

1  —  e 


(12) 


It  can  be  seen  from  Eq.  (11)  that  the  form  of  the  energy 
equation  for  constant  lag  two-phase  mixture  in  the  combustion 
chamber  is  the  same  as  that  for  pure  gas  phase  flow.  We  /£ 
refer  to  Cp  as  the  constant-pressure  specific  heat  of 
constant  lag  flow.  Note,  however,  that  this  is  not  the 
constant-pressure  specific  heat  of  the  two-phase  mixture,  but 
rather  a  symbol  introduced  so  as  to  be  able  to  write  the  equation 
in  the  form  of  that  for  pure  gas  phase.  The  total  temperature  TQ 
at  the  grain  head  is  equal  to  the  combustion  temperature  of  the 
grain. 


Similarly,  the  constant  lag  constant-volume  specific  heat 
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We  define 


£  —  (l~‘*)Cw+eCiL 
9  I -e+eK: 


Let 


i+rr*dL 

i  +*._£_6L 
1  — e 


(13) 


k  and  R  are,  respectively,  the  constant  lag  specific  heat 

ratio  and  gas  constant  of  the  combustion  chamber.  When 

K  =  L  =  0 ,  k  =  k,  R=  R  ,  and  we  have  pure  gas  phase  flow.  When 

9  t  , 

K  =  L  =  1 ,  1 _ 

t»k  *  - - - f  ^*(1* 

l+kr-2—t 

1-e 


corresponding  to  two-phase  equilibrium  flow.  Fig.  1  shows  the 
variation  of  k  and  R  with  particle  velocity  lag  (1-K) . 


Fig.  1.  Effect  of  particle  velocity  lag  on  k  and  R 


After  introducing  k  and  R,  we  obtain  from  Eq.  (11) 


(15) 
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In  the  above  equation, 


(16) 


X  is  termed  the  velocity  coefficient  of  the  constant  lag 
flow.  It  is  also  a  notation  adopted  for  the  purpose  of  obtaining 
a  solution  similar  in  form  to  that  for  pure  gas  phase,  and  does 
not  represent  the  actual  velocity  coefficient  of  the  gas  phase. 
acr  is  termed  the  constant  lag  critical  sound  velocity.  It 
varies  with  K  because  k  and  R  have  different  values  for  different 
values  of  K. 

P 

2.  Pressure  Ratio  pt 

___  If  fl  pV  p  £ 

Substituting  and  p#„#  “fTJ  into  Eq.  (5)  ,  we 

obtain 


(18) 


From  the  state  equation  of  the  ideal  gas, 


pt  P  T  • 
m~p~  *  r 

P*  P>  1  » 


/5 

(19) 

(20) 
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Substituting  Eqs.  (15)  and  (18)  into  the  above  equation,  we  get 

Pj L_,  1 

P#  1-faJ5  ^23 

P»  is  the  gas  phase  density  at  the  grain  head. 


4.  Total  Pressure  Ratio 


The  ratio  of  the  gas  phase  stagnation  pressure  Pg  at  any 
cross-section  of  the  grain  passage  to  the  pressure  P0  at  the 


grain  head  is 


On  the  other  hand, 


Ei.  Pjl  P 
p  ’K 


■(It)* 


Tg  is  the  gas  phase  stagnation  temperature  at  that  particular 
cross-section. 

The  gas  phase  parameters  at  the  same  cross-section  satisfy 
the  equation 

Cr'Tj+Aojjj—CrtTs 

It  follows  that 


T S  Ag  _ Uj _ u , 

29 

T  • 

Substituting  Eqs.  (15)  and  (16)  into  the  above  equation,  we 
obtain,  after  some  rearrangement. 


Tj_.  C, 
T'  ~^7; 
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Inserting  Eqs.  (18),  (b)  and  (c)  into  Eq.  (a),  we  have 

Pt  (l+an(l-|^ir)TS 


IT 

5.  Total  Temperature  Ratio 

Ts  Js.Ll 

T7  T7  T. 

Substituting  Eqs.  (15)  and  (c)  into  the  above,  we  obtain 


Therefore,  the  total  temperature  of  the  gas  phase  varies 
along  the  length  of  the  passage.  This  is  due  to  the  fact  that, 
in  the  flow,  the  particles  transfer  heat  to  the  gas  phase,  the 
gas  has  to  overcome  friction  due  to  the  particles  in  order  to 
do  work,  and  the  gas  phase  flow  process  is  not  an  adiabatic 
process. 

6.  Parameters  of  the  Particles 


The  flow  parameters  of  the  particles  are  calculated  from  the 
following  equations: 

vp~K  v, 

r,-(i-i)r.+Lr,  j 

n  e  j  (24) 

pf  ~(l-e)KPt 


Relation  between  J  and  \T 

XV  .  _ Jj 


From  conversation  of  mass, 


Hi  i,9*  fh  t 
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L-  A  £&k„  T  K 

1  — e  Af  P,  P*  ■  Zi**' 


(e) 


From  one-dimensional  steady  constant  two-phase  constant  lag 
nozzle  flow  theory^'  ^3] 


(f) 


In  the  above  equation 


r^i45^ 


C-l+S“  <*C*(l-iO+/n+(A-l)<5IZ)> 

r-vT^^r  _ 


Substituting  Eqs.  (e) ,  (f)  into  Eq.  (d) ,  and  noting 
that  PsL*P,a (It),  P.-p.gR,T„  ’  we  obtain 


'i* -i 


-  '  /  i  + 1  j 


/I 


(25) 


In  the  above  equation, 


,_L  r~i  I  a  ,r 

rvirrx^7pC 


It  can  thus  be  seen  that,  with  respect  to  constant  lag 
flows,  there  exists  a  definite  relation  between  the  XL  at  the 
tail  end  of  the  grain  and  the  throat-passage  ratio  J.  When  all 
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the  other  parameters  are  fixed,  \L  is  solely  determined  by  J. 
v  Equilibrium  Pressure  PQ  at  Grain  Head 


To  emphasize  the  study  of  the  effects  of  the  two-phase  flow 
in  the  combustion  chamber  on  the  performance  of  the  inner 
trajectory,  we  will  discuss  only  the  case  for  smaller  values  of  J 
(eg.  J  <  0.4).  In  this  case,  the  flow  velocity  is  relatively 
low,  and  the  effect  of  flow  is  not  serious,  so  that  one  may 
regard  the  combustion  rate  r  to  be  constant  along  the  entire 
length  of  the  passage,  and  equal  to  that  at  the  grain  head;  i.e., 
r  s  =  bP^4  .  It  follows  that  the  rate  of  increase  of  mass  of 
the  products  of  combustion  is 


ihi—  pr’bPi -At  (g 

From  the  assumption  of  instantaneous  equilibrium,  m^  =  mt. 

Inserting  Eqs.  (f)  and  (g)  into  the  above  equation,  we 
obtain  after  some  rearrangement 


In  the  above  equation. 


Mcv'wr'Jar^yj'-' 

*  _ L.  L  IT 

*  1-8  Ht 


1 

r 


r- 


2(A-1) 


(26) 
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VI .  Distribution  of  \  along  X 


Integrating  Eq.  (1)  from  0  to  x,  we  obtain 


As 


we  have 


Similarly, 


t*~f>Tr»Sx 

lh^A,iptV,  +pr  Vf)  —-^-yirptV, 

Y~A,p,v,=~pTr,Sx 
J^ZjAfP)  t«#  t  *  Prf  %S  I 


Dividing  the  last  equation  into  the  one  before,  we  have  /8_ 

PlV>  .  * -  » 

PtL«»L  1 

^-•e*  j£t.T 

_£i - - 

P*±-Il 

p» 

Inserting  Eq.  (21)  into  the  above  equation,  we  obtain  after  some 
rearranging. 


1+qfi 

It 


(27) 


It  can  be  seen  from  Eqs.  (15)  ,  (18)  ,  (21)  ,  (22)  ,  (25)  ,  (26) 
and  (27)  that  the  various  relations  pertaining  to  the  constant 
lag  two-phase  flow  in  the  combustion  chamber  are  similar  to  the 
corresponding  relationships  pertaining  to  the  pure  gas  phase. 
When  K  =  L  *  0 ,  we  obtain  the  relations  for  the  pure  gas  flow; 
when  K  *  L  =  1,  we  obtain  the  relations  for  the  two-phase 
equilibrium  flow.  Therefore,  the  pure  gas  phase  flow  and  the 
two-phase  equilibrium  flow  can  be  regarded  as  special  cases  on 
the  constant  lag  flow. 


VII .  Relationship  between  the  Lag  Constants  K  and  L 

From  the  mass  equation,  momentum  and  energy  equation  of  the 
particles  given  in  Ref.  [1] ,  we  can  derive 


dv_>  A(vt-v ,)  ...  S' 

dx  r. - ~W'T7 


Under  the  conditions  of  the  combustion  chamber,  the 

f  21 

particles  maintain  their  Stokes  flow  state  ,  and 


A-\  -JL>- 
2  r-rpmf  , 

p _ ^ 


From  the  definition  of  constant  lag  flow, 


dVp  =  Kdvg 
dTp  -  L*Tg 

Substituting  the  above  into  Eqs.  (28)  and  (29) ,  respectively,  we 
obtain 

dv.  A(l-K)  W,  S 
dx  “  7T5  X  A, 


From  Eq.  (11)  ,  we  have 


dT,  -  A«  dv, 
dx  ”  g&r*  dx  j 

r1,  > 

T*-T»  -  A°~2gS 7 
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Inserting  Eq.  (i)  into  Eq.  (h) ,  we  get 


/  9 


Comparing  the  two  expressions  for  ,  we  obtain 

ax 


A-iA**&-&)K 


(31) 


In  the  combustion  chamber  of  the  solid  propellant  rocket 

5 

engine,  A  and  B  are  on  the  order  of  magnitude  of  10 
2  gC, 


(1/sec)  , 


lVr^-«A 

('P  rrr  J £/ \  ..  & 

C,lV'A„  U  " K  )  gC, 


Therefore,  Eq.  (31)  can  be  simplified  as 


As 


we  have 


1 


1  —  K  r,  -4  , 

•2 gC,  -g-T  1 


~K~‘ 


B  2 


1 

gCpt 


1 


3^C^T  K  1 


(32) 


This  relationship  is  similar  to  that  between  the  lag 
constants  for  the  one-dimensional  two-phase  constant  lag  flow  in 
the  nozzle.  Thus,  the  chief  factor  affecting  the  relationship 
between  the  lag  constants  is  still  the  interaction  of  heat  and 
forces  between  the  particles  and  the  gas  phase.  The  effect  f  the 
particle  jet  velocity  Wp  can  be  neglected.  Eq.  (32)  shows  that 
when  and  Cpg  are  constant,  if  K  has  a  definite  value,  then  the 
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value  of  L  is  also  definite.  In  many  engineering  applications, 
f  21 

5  ,  and  0  <  K  <  1.  Hence  L  <  K,  or  1  -  L  >  1  -  K. 

This  shows  that  in  constant  lag  flows,  the  temperature  lag  of  the 
particles  is  greater  than  their  velocity  lag.  The  relationship 
between  L  and  K  is  as  shown  in  Fig.  2.  10 

0-8 
05 
04 
a  2 

0  02  04  04  0-d  10 

K 

Fig.  2.  Relationship  between  the  lag  constants  K  and  L 


VIII.  Effect  of  Particle  V elocity  Lag  on  Inner  Trajectory 
Performance 


In  reference  to  a  particular  solid  propellant  rocket  engine 
(e  *  0.26),  we  have  calculated  the  P-t  curves  and  the  flow  fields 
of  the  combustion  chamber  for  different  values  of  K. 

Fig.  3  shows  the  pressure-time  curves  for  different  values 
of  K.  It  can  be  seen  from  this  diagram  that  the  pressure  in  the 
combustion  chamber  with  velocity  lag  is  lower  than  that  of  the 
two-phase  equilibrium  flow.  Furthermore,  the  larger  the  particle 
velocity  lag,  the  lower  the  pressure.  The  reduction  in  pressure 
of  the  combustion  chamber  results  in  a  reduced  combustion  rare, 
causing  the  burning  time  of  the  grain  to  increase,  and  the  amount 
of  flow  to  decrease.  Table  1  gives  a  list  of  the  variation  with 
particle  velocity  lag,  at  t  =  0 ,  of  the  constant  lag  flow 
combustion  chamber  pressure,  grain  burning  time  and  the  ratio  of 
the  mass  flow  per  second  to  that  corresponding  to  the  two-phase 
equilibrium  flow. 

Fig.  4  shows  the  distribution  along  the  length  of  the  / 10 
passage  of  the  constant  lag  velocity  coefficient  A  and  the 
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Table  1 


l-K 

0 

0.1 

1.3 

0.3 

l  ^ 

D 

0.1 

0.7 

0.1 

0.9 

1. 000 

0.913 

0.916 

0.949 

0.131 

0.914 

fil 

1  O.S36 

■  1 

1.000 

1.003 

1.007 

'  i 

1.010  | 

1 

1.014  ^ 

1.011 

1.031 

' - i 

1.025  1 

1.030 

1 

1.035 

1.400 

0.994 

0.990  | 

0.917 

0.913 

5.  97*  j 

i 

i  1 

0.975  0.971  | 

0.567 

relative  gas  phase  parameters.  It  can  be  seen  from  the  diagram 
P  P  T  Ps 

that  -p-%  y~  and  ^  all  decrease  along  the  passage 

length,  while  X  increases  along  the  length.  The  trends  of 
variations  are  similar  to  those  for  the  two-phase  equilibrium 
flow.  The  values  are  slightly  lower  than  the  corresponding 
values  for  the  two-phase  equilibrium  flow,  but  the  differences 
are  very  small.  This  illustrates  that  the  particle  velocity  lag 
has  but  minor  effects  on  the  distribution  of  the  relative  values 
of  the  gas  phase  parameters  along  the  passage  length.  As  PQ 
decreases  with  the  increase  of  particle  velocity  lag,  the 
distribution  of  P,  Pg  and  *>*  along  the  passage  length  is  more 
severely  affected  by  the  velocity  lag.  It  can  be  seen  from  fig. 

5  that  the  values  of  these  parameters  are  all  smaller  than  those 
for  the  two-phase  equilibrium  flow,  and  the  greater  the  velocity 
lag,  the  lower  these  values.  However,  the  particle  velocity  lag 
has  very  little  effect  on  the  distribution  of  T^. 

The  distribution  of  the  gas  phase  velocity  v^  and  the 
particle  velocity  vp  along  the  passage  length  under  different  lag 
conditions  at  t  =  0  is  shown  in  Pig.  6.  All  the  values  increase 
along  the  passage  length,  characteristic  of  the  introduction  of 
mass.  However,  K  has  opposite  effects  on  v^  and  vp.  It  can  be 
seen  from  the  diagram  that  in  the  presence  of  a  lag,  vp  has  a 
lower  value  than  that  of  the  two-phase  equilibrium  flow:  the 
greater  the  lag,  the  smaller  vp.  On  the  other  hand,  vg  has  a 
higher  value  than  that  for  the  two-phase  equilibrium  flow;  the 
greater  the  lag,  the  greater  v  .  This  is  because  when  (1  -  K) 
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increases,  both  k  and  R  increase,  causing  a  to  increase,  while 

the  particle  velocity  lag  has  a  very  small  effect  on  X. 

Therefore,  when  (1  -  K)  increases,  v  increases.  From  the  stand- 
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point  of  physical  significance,  the  value  of  K  actually  reflects 
the  size  of  the  particles.  A  large  velocity  lag  (small  K) 
signifies  a  large  size  of  the  particles.  We  know  from  Eq.  (30) 
that  Ao c-Lf  .  Hence,  for  a  large  lag,  A  is  small,  the  drag 

X  =  A(Vg  -  Vp)  exerted  on  the  gas  phase  by  a  unit  mass  of 
particle  is  small,  and  as  a  result,  the  gas  phase  velocity  is 
large.  Of  course,  under  such  conditions,  the  force  with  which 
the  gas  phase  carries  the  particles  along  is  also  small.  The 
particles  have  a  low  acceleration,  and  their  velocity  is  low. 

The  distribution  of  the  temperature  and  density  of  the 
particles  along  the  length  of  the  passage  is  shown  in  Fig.  7. 

The  decrease  of  *>/*,  and  Tp  along  the  passage  is  very  small. 

The  particle  velocity  lag  has  an  extremely  small  effect  on  the 
distribution  of  Tp,  but  a  rather  large  effect  on  the  distribution 
of  Pr  .  It  can  be  seen  from  the  diagram  that  the  value  of  pf 
in  the  presence  of  lag  is  higher  than  that  for  the  two-phase 
equilibrium  flow,  and  the  greater  the  lag,  the  greater  the  value 
of  . 


The  results  of  computations  show  that  the  total  temperature 
Tg  of  the  gas  phase  changes  very  slightly  along  the  passage.  For 
instance,  for  1  -  K  =  0.9  and  t  =  0,  tQ  =  3042. 4 °K  at  the  head, 
while  Tc  =  3042. 46 °K  at  the  end  of  the  arain  passage.  The 

OjL  '  £) 

increase  is  only  0.06°K.  Actually,  as  in  Eq.  (23)  is  very 

close  to  unity,  Tg  can  be  considered  as  being  constant  along  the 
passage  in  the  computation.  As  a  result,  Eqs.  (22)  and  (25)  can 
be  simplified  to 


r»‘";  7 


(22') 


(25*  ) 
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IX.  Conclusion 


/II 


We  can  conclude  from  the  above  discussion  that: 

1.  The  fundamental  relations  for  the  one-dimensional  two- 
phase  constant  lag  flow  in  solid  propellant  rocket  engines  are 
similar  to  those  for  the  pure  gas  phase  flow.  Pure  gas  phase 
flow  and  two-phase  equilibrium  flow  can  both  be  regarded  as 
special  cases  of  constant  lag  flow.  Therefore,  the  method 
presented  in  this  paper  for  inner  trajectory  computations  has 
greater  generality. 

2.  Owing  to  the  effect  of  the  particle  velocity  lag,  the 
pressure  in  the  combustion  chamber  decreases,  the  burning  time  of 
the  grain  is  lengthened,  and  the  nozzle  flow  is  reduced. 

3.  Two-phase  flow  has  a  relatively  large  effect  on  the  flow 
field  of  the  combustion  chamber.  With  increasing  particle 
velocity  lag,  the  gas  phase  flow  rate  increases,  pressure,  total 
pressure  and  gas  phase  density  decrease,  particle  velocity 
decreases,  and  particle  density  increases.  However,  the  particle 
velocity  lag  has  little  effect  on  the  temperature  of  the  gas 
phase  and  the  particles. 


Fig.  3.  Pressure-time  curves  Fig.  4.  Distribution  of 

relative  values  of  gas  phase 
parameters  in  two-phase  flow 
Key:  1)  two-phase  equilibrium  flow 
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Fig.  5.  Distribution  of  gas  phase  parameters  in  two-phase  flow 


Fig.  6 
Key: 


Fig.  7 


.  Distribution  of  velocity  of  the  gas  and  the  particles 
two-phase  flow 

1)  two-phase  equilibrium  flow 


.  Distribution  of  particle  parameters  in  two-phase  flow 
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